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Abstract. We study the three dimensional Ginzburg-Landau model of su- 
perconductivity. Several 'natural' definitions of the (third) critical field, Hc 3 , 
governing the transition from the superconducting state to the normal state, 
are considered. We analyze the relation between these fields and give con- 
ditions as to when they coincide. An interesting part of the analysis is the 
study of the monotonicity of the ground state energy of the Laplacian, with 
constant magnetic field and with Neumann (magnetic) boundary condition, in 
a domain f2. It is proved that the ground state energy is a strictly increasing 
function of the field strength for sufficiently large fields. As a consequence of 
our analysis we give an affirmative answer to a conjecture by Pan. 



1. Introduction 

In the whole paper £1 C R 3 will be a bounded simply connected domain 
whose boundary is connected and Lipschitz continuous. 

Let j3 be the constant magnetic field along the z axis: (3 = (0, 0, 1). The Ginzburg- 
Landau functional in three space dimensions is given by 

£[i>,A] =E K , H [iP,A] =^{b K HAVf -« 2 M 2 + yM 4 }^ 

+ k 2 H 2 I \cmlA- /3\ 2 dx , (1.1) 

JUL 3 

with tjj G W 1,2 (£l; C), A in the space Hp div that we will define below, and where 
Pa = {—iV + A). Notice that the second integral in (jl.ip is over the entire space, 
M 3 , whereas the first integral is only over the domain f2. 
Let F be the vector potential 

F(xi,x 2 ,x 3 ) = ±(-x 2 ,xi,0) . (1.2) 

Formally the functional is gauge invariant. In order to fix the gauge, we will 
impose that vector fields A have vanishing divergence. Therefore, a good choice 
for the variational space for A is 

H FAiv = F + H div i (1-3) 

where 

Hl iv ={Ae H 1 ^ 3 ,^ 3 ) | div A = 0} . 

We use the notation ij 1 (R 3 ) for the homogeneous Sobolev spaces, i.e. the closure 
of C5°(IR 3 ) under the norm 

/ -> 11/11^ = IIV/Hl* . 
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We will recall below the fact that any square integrable magnetic field B, i.e. any 
vector field B £ i 2 (R 3 ;R 3 ) with divB = in the sense of distributions, can be 
represented by a vector field A e H^ iv . 

Minimizers, (ip, A) <E W 1,2 (VL) x Hp div , of the functional £ have to satisfy the 
Euler-Lagrange equations: 

pLfA^ = « 2 (l-M 2 )^ in n, (1.4a) 
curl 2 A = { - ^h^pVip - ^V^) + \iP\ 2 A}l n (x) in W 3 , (1.4b) 

(p K HAlb)-N = on 90, (1.4c) 

where N(x) is the unit interior normal at the boundary. 

It is not completely standard (but see [GiJ for the analysis of this case) but 
rather easy to prove that, for all n, H > 0, the functional £ Kj h has a minimizer. 
It is a result of Giorgi and Phillips, [GiPhj . that for k fixed and H sufficiently 
large (depending on k), the unique solution of ()1.4|) (up to change of gauge) is the 
pair (V')A) = (0,F). Since ip is a measure of the superconducting properties of 
the state of the material and A is the corresponding configuration of the magnetic 
vector potential, the result of Giorgi and Phillips reflects the experimental fact that 
superconductivity is destroyed in a strong external magnetic field. 

We define the lower critical field, H_ C[J as the value of H where this transition 
takes place: 

Hjq (k) :— inf{_ff > : (0,F) is a minimizer of £ k ,h} ■ (1.5) 

However, it is far from obvious from the functional that the transition takes place 
at a unique value of H — there could be an interval of transitions back and forth 
before the material settles definitely for the normal state, (0, F). Therefore, we 
introduce a corresponding upper critical field 

Hc 3 (t) := m£{H > : for all H' > H, 

(0, F) is the unique minimizer of £ k ^h'} ■ (1-6) 

Notice that our space H F div fixes the choice of gauge. 

It is the objective of this paper to give conditions under which these two defini- 
tions of the critical field coincide (the transition being then a sharp phase transition 
at a precise value). More precisely, we will study the relation of the fields defined 
above to analogous local ones given purely in terms of spectral data. 

First let us recall the result on the asymptotics of the critical field. In |LuPa4[ 
IPan2] ) it was proved that there exists a universal constant Go ~ 0.59 — defined in 
(|3.1[) below — such that for all (smooth, bounded, simply connected) Slcl 3 there 
exists a constant C > such that for k large 

\H C3 (k)-£- I <C^ 2 . (1.7) 

Here Hc 3 {k) denotes either Hc 3 {k) or H_q (k). 

The local fields are determined by the values where the normal solution (0, F) is 
a not unstable local minimum of £ Kj Ht i- e - 

H 1 cI(k) := M{H > : for all H 1 > H, Hess£ K>H '\, F) > 0} , 

H 1 S c 3 (k) := mi{H > : Ress£ K>H | (0F) > 0} . (1.8) 
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Since the Hessian, Hess£ K ,#, at the normal solution defines the quadratic form 

(<£,a)i-> / + nHF)<f>\ 2 - K 2 \cj>\ 2 dx + (kH) 2 [ \cmla\ 2 dx, (1.9) 

Jn Jr 3 

we get the equivalent definitions given by 

H 1 c 3 (k) = vaf{H > : for all H' > H, X^kH') > k 2 } , 

H}£ c 3 (k) = wf{H > : Xi(kH) > k 2 } . (1.10) 

Here Xi(B) is the the lowest eigenvalue of the magnetic Neumann Laplacian H(B), 
i.e. of the self-adjoint operator (with Neumann boundary conditions) associated to 
the quadratic form 



W 1 > 2 {n)3 Q B {u) := [ \p BF u\ 2 dx 

Jn 



(1.11) 



In other words, H-(B) is the differential operator p 2 BF with domain {u £ W 2 ' 2 (il) : 

N ■ PBFu\on — 0}. The operator Tt(B) clearly has compact resolvent. 

The analysis in |LuPa4[ IPan2| implies that f| 1 . T|) remains true for the local fields, 

i.e. 

\H^(k)-^-\<CkV\ (1.12) 

where Hq^(k) denotes either (k) or Hj^ (k). 

Before stating our first main result let us give a third possible definition of 
the critical fieldaj. In two dimensions the magnetic field energy is usually given 
as an integral over the domain only. By analogy, one finds the following slight 
modification of the GL-functional 

£ 1Jlod [tp, A] = £™ H d {^ A] - J^{\ PkH a^\ 2 - « 2 H 2 + yH 4 } dx 

+ k 2 H 2 I |curlA-^| 2 dx . (1.13) 
Jn 

Here (i/>, A) £ W 1 ' 2 (ri;C)xW 1 ' 2 (Q;W. 3 ). Using the gauge invariance of the problem, 
we can and will assume A to be restricted to the smaller space 

Hl iv {Q) := {A e VF 1,2 (n;R 3 ) | div A = in fi , N ■ A = on dfl } . (1.14) 

This modified functional leads to a new set of possible values for the critical field. 
Let Fq denote the vector potential in H^ iv (O) generating the constant magnetic 
field in f2, i.e. 

divFn =0 

curlFo =(i 

The new fields are given analogously to the previous ones. 



inQ, N-F n = on dQ . (1.15) 



Hc° d (n) := mf{H > : for all H' > H, 



(0, Fq) is the unique minimizer of S^jji} , 



H^f{n) := mi{H > : (0, F n ) is a minimizer of S^} . (1.16) 



1 Our £ is probably the physically most correct of the two. In the literature both £ and £ mod 
are considered. For example, [L"uPa4 study £ mod , but in the follow-up paper Pan2 , it is £ that 
is taken as the definition of the functional. 
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The leading order asymptotics (|1.7p also holds for i?™° d and H^° d (n). 

We also state the Euler-Lagrange equations for stationary points of the modified 
functional 

cur^A^-^^W- W) + M 2 A J ' 1 U&) 



on dn. (1.17b) 



(PkjtaVO • N = 
(curl A - (3) x AT = 
Remark 1.1. 

FFe do not have to define the local fields corresponding to £ mod — these coincide 
with the previously defined local fields. To see this, notice that since Q is simply 
connected and curlF|^ = curlFn = [3 there exists a gauge 4> on Q such that 

P|n = Pfi+V0. 

Therefore, the magnetic Neumann operators (— iV + -BF) 2 and (— iV + BFq) 2 on 
L 2 (Q) are unitarily equivalent. In particular, 

inf Spec(-iV + BF) 2 = inf Spec(-iV + BF n ) 2 ■ 

Our first main result (combining Proposition 11.21 and Theorem 1 1 . 3|) is that all 

the critical fields above are contained in the interval [H 1 °,^(k), Hq 3 (k)], when k is 
large. We first observe the following general inequalities. 

Proposition 1.2. 

The following general relations hold between the different definitions of Hc 3 ■' 

H^ c M<Ec 3 (K)<H^ 3 d ( K ), (1.18) 

HcI(k)<HcM<HcT(k)- (1-19) 
For large values of n, we have a converse to Proposition 1 1.21 
Theorem 1.3. 

There exists kq > such that for k > kq, 

H% c 3 (K)=H jC3 (K)=m° d ( K ), (1.20) 
H } ° c U)=H C3 (K)=Hc° 3 d (K). (1.21) 



An important consequence of Theorem 1 1.31 is that in order to obtain an asymp- 
totic expansion of Hc 3 (k) for large values of k, whichever the definition, one only 
has to consider the linear problem of determining Xi(B). 

We state the following useful corollary to Theorem 1 1.31 



Corollary 1.4. 

Suppose that there exists Bq > such that B t— > Xi(B) is strictly increasing for 
B > Bq. Then there exists k > such that for k S [k 0i oo) one has 

= K C3 (n) = H™ d ( K ) = H V °1{k) = HcM = i?c 3 ° d («) . (1.22) 



Remark 1.5. 

The monotonicity of Ai (B) can be proven if one has a sufficiently good asymp- 
totic expansion for large B. For domains with smooth boundary (as considered 
here), we have, Xi(B) = QqB + o(B), to leading order. It follows from |FoHe2| 
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Proposition 2.2] (the proof is independent of dimension) that if one can prove an 
asymptotics up to order o{\), i.e. 

M 

Ai (B) = ® B + ^2 ajB 7i + o(l) , (1.23) 

with M G N, aij G ffi, 7j € [0, 1), t/ien Ai(-B) is monotonically increasing for large B. 
However, to the authors ' knowledge, there is no example in three dimensions, where 
an asymptotics as precise as (| 1 .23(1 is known (except for cylindrical, i.e. effectively 
two-dimensional cases). Nevertheless, as we will see below, one can modify the proof 
of |FoHe2| Proposition 2.2] in order to get the monotonicity of Ai as a consequence 
of a less demanding asymptotics combined with localization estimates on the ground 
state eigen] 'unction. 

The analysis of Hc 3 described above clearly calls for a clarification of the mono- 
tonicity of the function B i— > Ai(£>). In the two-dimensional situation this has 
been analyzed in [FoHe3|, IFoHe4j . We will here carry out the similar study of the 
three-dimensional case. We will prove (under a generic assumption on the domain 
f2) that the mapping B i— > \i(B) is monotonically increasing for sufficiently large 
values of B. 

We will work under the following geometric assumptions. 
Assumption 1.6. 

The set of boundary points where [3 is tangent to d£l, i.e. 

T:={xe dn\(3-N(x) = 0}, (1.24) 

is a regular submanifold of dfl. 

Therefore, T is a disjoint union of regular curves. We choose an orientation on 
each such curve, and define the normal curvature at the point x G T by 

k n (x) :=K x (T(x)AN(x),f3). (1.25) 

Here K denotes the second fundamental form on dfl, and T(x) is the oriented, unit 
tangent vector to T at the point x. We assume that 

Assumption 1.7. 

k n ( x ) ^ 0, Vx G T. (1.26) 

and that 
Assumption 1.8. 

The set of points where (3 is tangent to T is finite. 

Assumptions 11.61 11.71 and 11.81 are clearly generically satisfied. They are for in- 
stance satisfied for ellipsoids, whereas a domain containing a cylindrical boundary 
piece with axis parallel to [3 will violate these assumptions. 

We will need the known two-term asymptotics of the ground state energy of 
H(B). The following result was proved in HeMo5J (the corresponding upper bound 
was also given in Pan2j and a less general geometric situation was studied in 
[HeMo4] ). 
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Theorem 1.9. 

There exist universal constants 00,^0,00 > (to be defined in (|3.1|) . (|3.2p . (|3.6|) 
below) such that if il satisfies Assumvtions \77Sl |J.7| awrf li.ffl £/ien 

Ai(S) = e S + 7oBi (1.27) 

for some r] > 0. 

Here 70 is defined by 

70 := inf 70(1), (1.28) 

7o(x) := 2- 2 ^ d 1 /3 \k n (x)\y^6 + (1 - S )\T(x) ■ ft 2 ) 1 '* . (1.29) 

Notice that o"o < 1, so there is no problem with the third root in (|1.29|) . Com- 
bining this result with Corollary 1 1.41 one gets a two-term asymptotics for Hc 3 {k): 

Corollary 1.10. 

Suppose that ft C M 3 satisfies A ssumvtions [7751 11.7| anrf li.ffi f/ien one /inds 

HcM (£- - %®o 2/3 k 1/3 ) = o(n^) , (1.30) 

where Hc 3 (k) denotes any of the six different (upper or lower) critical fields defined 
above. 

Corollary 11.101 is an affirmative answer to a conjecture in |Pan2] . however, the 
conjecture is stated without the geometric Assumptions 1 1 . 6l 1 1 . 71 and 11.81 
Proof of Corollary \1.1(A 

By Proposition 11.21 and Theorem 11.31 it suffices to prove that H% c 3 (k) and fl«£(«0 
have the asymptotics given by p.30p . But this follows easily from (|1.27p . □ 

Finally we prove that that X\(B) is increasing for generic f2, which implies that 
the different critical fields coincide, cf. Corollarv ll.4l 

Theorem 1.11. 

Let fl C R 3 satisfy A ssumvtions \1.6l |J.7| and \1.8[ Let {T±, . . . , r„} be the collection 
of disjoint smooth curves making up T. We assume in addition that : 

For all j there exists x S Fj such that jo(x) > 70 ■ (1-31) 
Then the directional derivatives A' x ± :— lim^o± > sxist and satisfy 

lim a; + {b) = lim a; _{b) = e . (1.32) 

B — >oo ' B — >oo ' 

In particular B i— ► Xi{B) is strictly increasing for B sufficiently large. 



2. The analysis of Hc 3 

In this section we give the proof of Theorem ll.3l We aim at giving a simple proof 
without refering to long technical papers. Therefore, some of the statements in this 
section will not be the best possible ones known in the litterature. In particular, we 
avoid the use of the elliptic estimates obtained by 'blow-up' techniques (cf. |Pan2j 
and see also |LuPal|. [Mm\ IFoHe4j ). 
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2.1. Basic estimates. 

In this section we mainly recall a number of results that will be useful in the proof 
of Theorem 11.31 First we need some weak control of Xi(B). 

Theorem 2.1 (First order eigenvalue bound). 

Let Slcl 3 be bounded and simply connected with smooth boundary. Then 

Ai(B) = Q a B + o(B) . 

Much more precise estimates are proved in }LuPa4i [HcMo5 , Pan2] under addi- 
tional assumptions (cf. Theorem 1 1.9[) . but we want to stress that, for the proof of 
Theorem II .3) . we only need very weak spectral information. 

It is a standard consequence (cf. |DGPj ) of the maximum principle that a mini- 
mize iQ (i/>, A) of £ K . H or satisfies 

< 1 . (2.1) 

The normalization of our functional £ Ki h is such that £ Kj #[0, F] = 0. So any 
minimizer (tp,A) will have non-positive energy. Therefore, the only negative term, 
— ^H^Hii m the functional has to control each of the positive terms. This leads to 
the following basic inequalities for minimizers, 

WPKHAlPh < «||V||2, (2-2) 

F||curlA-/3|| 2 < ||V|| 2 . (2.3) 

Furthermore, using (|2.ip . 

MI < Mia. (2-4) 

The same inequalities remain true for minimizers of • 
Finally, we need elliptic estimates for the curl-div system: 

Theorem 2.2 (Ellipticity of the curl-div system). 

There exists a constant C > such that for all (magnetic fields) b £ L 2 (R 3 ,R 3 ) 
with divb = 0, there exists a unique a £ 7J 1 (R 3 ,K 3 ) such that 

curl a = b , div a = . 

This solution satisfies the estimate 

Hallux < C \\h\\ L * . (2.5) 

Proof. 

An argument for this standard result is given in [GiPh . It is based on the elemen- 
tary fact that, for / e Cg° (R 3 ;R 3 ) one has 

WfWfr = I |div/| 2 + |curl/| 2 dx. (2.6) 

With r(x) = 47^1 being the fundamental solution of the Laplacian, the desired 
solution is (formally) a = —curl (r * b). □ 

Proposition 2.3. 

Let 2 < p < 6 and let fl C M 3 have bounded measure. Then there exists a constant 
C p > such that for all b 6 L 2 (M 3 ,R 3 ) with divb = the solution a given in 
Theorem \2.2\ satisfies the estimate 

IHlppi) < C p \\h\\ L 2 (m) . (2.7) 



2 More generally, the inequality 1 12. Il l is valid for all stationary points. 
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Proof. 

By (|2.5|) and the standard three dimensional Sobolev estimate 

\\f\\L* m < C Soh ||/||^ , V/ e H\M. 3 ) , (2.8) 

the desired estimate holds forp = 6. Since SI has finite measure, Holder's inequality 
implies that HaH^n) < C||a|| L 6 (n) , for p < 6. □ 

The estimates we actually use below are the ones corresponding to Theorem 12. 21 
and Proposition 12.31 for vector potentials a and magnetic fields b defined only on 
51. Recall the definition of H^ iv (SI) in (|1.14[) . The next theorem is standard. It is 
a summary of the results (in the simply connected case) of the discussion in |Te[ 
Appendix 1]. 

Theorem 2.4. 

Let SI C M 3 be smooth, bounded and simply connected and having a finite number 
of connected boundary components denoted by S\, . . . , S n . Then the space of L 2 - 
magnetic fields, i.e. the image o/W /1,2 (0) under the operator curl , satisfies 

S:=ciMW ll2 (Q) = {be£ 2 (fi)|divb = 0, / b • TV dS = , Vi } . (2.9) 

J Si 

In particular, if n = I, we have 

B = {be L 2 (Sl)|divb = 0} . (2.10) 

Furthermore, for all h £ B, there exists a unique a £ -ff^iv 0^) satisfying 

curla = b. (2.11) 

Finally, there exists a constant C > ( depending only on 51 J such that for all 
b £ B, the solution a £ H\ iv (51) to (|2.1ip satisfies 

||a|| w i,2 ( n) < C \\h\\ L * {n) . (2.12) 

We get an analogous estimate to the one in Proposition 12.31 with unchanged 
proof. 

Proposition 2.5. 

Let 2 < p < 6 and let il C M 3 be bounded and simply connected with smooth 
connected boundary. Then there exists a constant C p > such that for all b £ 
L 2 (f!),R 3 ) with divb = the solution a given in Theorem \2.4\ satisfies the estimate 

l|a|U*(Q) < C p \\b\\ L 2 {n) . (2.13) 



2.2. Weak decay estimates. 

In this subsection we prove the following decay estimate in the variable normal to 
the boundary. Much more precise (Agmon type) decay estimates exist but they 
depend on a much longer analysis of the Ginzburg-Landau equations (see |Pan2j ). 
The estimate below is very robust and sufficient for our purpose. It extends [ BonFo] 
who treat the two-dimensional case. 

Theorem 2.6 (Weak normal decay estimate). 

Let n C M 3 be a bounded domain with Lipschitz boundary. Then there exist two 
positive constants C andC, such that, if (tp, A) K> # is a minimizer of £ K} h or , 
with 



k(H-k) > 1/2, 



(2.14) 
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then 

Hi<cf Wx)\ 2 dx< - C (2.15) 

J{ x /k(H-k) dist(x,90)<l} ^/k{H — k) 

and 

H\\mn)<C p [K{H-K)Y P -^U\\ LP{n) . (2.16) 

Proof. 

There is no real modification for the 3D case in comparison with p BonFbj but 
we repeat the proof for completeness. The last inequality in (|2.15[) is an easy 
consequence of (|2.1j) . since there exists a constant C\ > (depending only on f2) 
such that 

meas{a; : dist(a;,9fi) < A} < C x \ , VA e (0,2] . (2.17) 
Let x G C°° 0^) be a standard non-decreasing cut-off function, 
X = l on [l,oo), X = on (—oo,l/2). 



Define A := 1/^/k{H - k) and xa : & -> K by 

Xa(«) := x(dist(x,50)/A). 

Then xa is a Lipschitz function and suppxA C fi. Combining the standard local- 
ization formula and (|1.4a[) . we find 

\p K HA(xxn 2 dx - f IVxaHV'I 2 dx = ^(xl^,n KH A^) 



= k 2 J \ X x4>\ 2 dx - k 2 J xiWdx. 

(2.18) 

Note also that, by integration by parts, one has, since xx^P has compact support, 
the following important though elementary inequality, 

/ \p K HA(xxn 2 dx > kH f (curlA^XAVf (2.19) 
Jn Jn 

Then, using (|2.3|) . we get 



/ 

J si 



P k ha{xx4>)\ 2 dx > K H\\x^\\l - nH\\cm\A - /?| 1 2 1 1 Xa^ II 2 

> kH\\xxiI>\\1 ~ \U\\l ~ « 2 HXA^Ilt (2.20) 



Using ([23]) and ((231), we get from (|2~Tg|) and ([2~2"TJ|) that 

K(H-K)||XA^Il2 

< \U\\1 + Hx'llLA- 2 / \i>(x)\ 2 dx + K 2 !{x\-xl)Wdx. 

Notice that the last integral is negative and we thus find by splitting HV'Hl 

{/ S (fT-K)-l/4}||xAV'||2 

< (H^'llSoA- 2 -h 1/4) / \^{x)\ 2 dx. 

By assumption 

k(H - k) - 1/4 > k(H -k)/2. 
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Moreover the conditions on x an d k(H — n) imply that 

||x'||LA- 2 + l/4<(||x'llL + l)A- 2 . 

Thus, 

\M\\l<2(\\x'\\l + l) / mx)\ 2 dx. (2.21) 

J{dist(x,an)<\} 

Consequently, 

llV'll^<4(||x'l| 2 ; + l) / \m\ 2 dx. (2.22) 

./{dist(x,9fi)<A} 

This finishes the proof of (|2.15|) . 

The L p estimate f|2 . 16f) is a consequence of the first inequality in (|2.15j) and 
Holder's inequality. □ 



2.3. Equal fields. 

We first give the easy proof of Proposition 1 1.21 

Proof of Proposition 

The inequality Hq^(k) < H c (k): 

Suppose H < H^(k). Then Xi(kH) < k 2 . Let ip be a ground state for H(kH). 
We use, for rj > 0, (r)ip, F) as a trial state in £ K: h, 

F] = (Xi(nH) - n 2 ) v 2 U\\h m + Y^UWUn) ■ 

Since Xi(nH) — k 2 < 0, we get £ K} h['<1' i P; F] < for n sufficiently small. Thus (0, F) 
is not a minimizer for £ Kj h- Since H < H_§^(k) was arbitrary, this proves that 

The inequality H c? (n) < H]S,° d (K): 

Suppose H < H_c 3 ( K )- Then there exists (V>,A) € W 1 ' 2 (fl) x div such that 
£ Kt jj[^j,A] < 0. By restriction A defines an element A 6 W 1,2 (f2,R 3 ) and we get 
the following simple inequalities, 

£™ H d [^,A] <£ K , H [>I>,A] <0. (2.23) 

Thus (0, Fn) is not a minimizer for £™#. Since H < H c (k) was arbitrary, this 
proves that H C3 (k) < H^° a {k). This finishes the proof of ([Tig]) . 
The inequality H C3 (k) < Hc 3 (n): 

Suppose H > Hc 3 (k). Then (0, F) is the only minimizer of £ Kj h- In particular, for 
all s £ R and all cj>, A, 

£ K)ff [#,F + sA] >£ K , H [0,F] = 0. 
This implies that Hess£ K .# |, Q p ^ > 0. Since H > Hc 3 (n) was arbitrary this proves 
that Hclin) < Hc 3 (k). 
The inequality Hc 3 (k) < i/™° d («;): 

Let H > H C3 («) and let ft£>, A) € W 1,2 (n) x i? J div be a minimizer of £ k ,h- Then, 
since £ Ki ir[0, F] = 0, we find, with A being the restriction of A to f2, 

0<£™ d [^i]<£ K ^,A]<0. 
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Furthermore, since H > i?™° d («;), the first inequality can only be an equality if 
i/j = 0. But then we find that 

= £ Kt jj[ip, A] = (kH) 2 / |curlA-/3| 2 da;, 

i.e. curl A = f3, which implies that A is gauge equivalent to F. Thus, for all 
H > ( K ) ^ ne om y minimizer of £ k ,h is the normal state (0,F). Therefore 

we have proved the inequality Hc 3 (k) < ff™ d (K). This finishes the proof of 

unsj). □ 

The proof of (|1.21| is essentially identical to the proof of the corresponding 
statement in two dimensions given in |FoHe2j . The small difference occurs in the 
use of Proposition I2.5I We give the proof here for completeness and to emphasize 
that it is remarkable that this argument is sufficiently robust to be generalized to the 
three dimensional setting where the spectral asymptotics is much more complicated 
than for planar regions. 

Define the sets 

Af(n) := {H S R | £ k .h has a non-trivial minimizer} , (2.24) 
7V mod (K) := {H e R | £™°h has a non-trivial minimizer} , (2.25) 
Mi oc (k) -~ {H eR\\i(nH) < k 2 } . (2.26) 

By evaluating the functional £ k ^h (resp ) in the state (77^, F) (resp (rjip, Fn)), 
with rj small and tp being the ground state ofH.(nH), one gets the inclusions, 

Moc(k) C AT(k) , Moc(k) Q Af mod (K) . (2.27) 



This, of course, is analogous to Proposition 11.21 We will prove that the converse 
inclusion holds for large k, so we get 

Theorem 2.7. 

Let n be bounded and simply- connected with smooth boundary. Then there exists 
kq > such that for all k > kq we have 

Moc(k) = 7V(k) = J\f mod (n) . (2.28) 

Of course, p.27p implies that we only need to prove inclusions in one direction 
in (|2.28p . We will only prove A/" mod (ft) C Ni oc (k), the similar inclusion for N(n) 
being proven in exactly the same way. 

Lemma 2.8. 

Let c > 0. Then there exists k > such that if 

H-k>ck, (2.29) 
and (tp, A) K> _ff is a nontrivial minimizer of £™ 1( }f with k > Kq, then 

ft 2 — Xi(kH) > . (2.30) 
Actually, we have the more precise estimate 

< ^\\ni-Q. H .m < (1 + 0{Krkw _ Xi{KH)) {2M) 
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Proof. 

Recall the rough asymptotic estimate 

\H^ d (n)--^\<CV^, (2.32) 

stated in the introduction for any definition of Hc 3 ■ 

By definition of i?™° d (K), nontrivial minimizers only exist below _ff™° d (ft), so 
we may assume that 

(l + c) K <H<Hc° d (K) . 
Since (ip,A) is non-trivial, we get that : 

K 2 H\\ 2 2 > Q k ha.VI>] • (2.33) 

We define 

A:= k, 2 \\tP\\ 2 2 -Q kHA [tP] . (2.34) 
Notice, that the GL-equation gives 

U\\i = 4' ( 2 - 35 ) 

hi 

We are in a situation where Theorem 12.61 can be applied. Therefore, we get by 
([2~T6) with p = 4 

Wh<CK-i\\i>\U. (2.36) 

Coming back to (|2.35p . we get 

\\i/j\\ 2 < Cn-iAi . (2.37) 

We now replace A by Fq in the quadratic form. Denote a = A — ¥q and b = 
curl (A — Fq). Then the Cauchy-Schwarz inequality implies that 

< A < [hi 2 - (1 - P )X 1 (kH)] HV-lla + P^\nH) 2 [ \aifj\ 2 dx , (2.38) 

for all < p. 
Proposition 12.51 implies 

||a|U*(Q) < C\\b\\ L * {Q) , 

Therefore, 

(kH) 2 \\&\\ 2 < C(hiH) 2 \\cmlA - < CA . (2.39) 

Here we used that £™]f[ip,A} < to get the last estimate. 
We now insert (j2~39j) . (j2~35]) . and (|2~37|) in (|2~38|) . 

< A < [k 2 - (1 - p)\i(kH)] Ml + P -\kH) 2 \\*\\ 2 U\\ 2 



< [k 2 - Xi{kH)] \\tp\\l + Cp\i(KH)A^K-i + Cp- l A— (2.40) 

Upon choosing p = VAk"^, and using that Xi(kH) < Ck 2 , we find, for another 
constant C, 

< A < [k 2 - Xi(kH)] \\ip\\l + CAk~^ . (2.41) 
When n is so big that CkT^ < 1, we therefore get 

< (1 - CV^A < [k 2 - A x («;#)] . (2.42) 
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Since tp cannot vanish identically for a non-trivial minimizer, this shows both (pOO]) 
and fOTj) . □ 

3. MONOTONICITY OF Al 

3.1. Spectral theory. 

In this subsection we will recall a few results of the spectral theory of some impor- 
tant model operators. This will permit to define precisely the constants 0o, S and 
£>o appearing in the main statements. 

3.1.1. The model on M 2 ^. 

Consider first the self-adjoint operator P R 2_ defined on L 2 (R+) (with R+ := {(s, t) S 
R 2 : t > 0}) by the quadratic form, 

{uel 2 (l+) : (-iV+P E 2)tt6L 2 (ll)}9uw / l(-iV + F R 2)w| 2 dsdt. 

Jm 2 + 

Here F R2 (s,t) := |(-t,s). 

The first constant that we need to define is 

6 := inf Spec P R 2 . (3.1) 

It is known that | < Go < 1- 

A unitary transform reduces the study of P R 2_ to the study of the Neumann real- 
ization of D\ + (t - s) 2 on L 2 (R+), which again reduces to the study of the family 
of operators (s 6 R), 

f)( S ) :=D 2 + (t-s) 2 , 

on L 2 (R_|_) with Neumann boundary conditions. In particular, with the notation 
fi(s) :~ inf Spec f)(s), we have 0o = inf sS R /i(s). One can prove that this infimum 
is a unique, non-degenerate minimum. I.e. there exists a unique sq S R such that 
Bo = m( s o), and 

*o := ^ > 0. (3.2) 

Actually, 8 € (0, 1), and an alternative expression for it was derived in [FoHel] . 
For more information on the spectral theory of P R ^ see [DaHel IHeMo2j . 

3.1.2. The model on R\. 

The next model operator to study is the case SI = R^_ := {(xi, X2, X3) el 3 : x\ > 
0}, (3 = (sm9, cosd, 0), B = 1 in p. lip . In other words, and with a particular 
gauge choice, it is the Neumann realization of the differential operator 

P(Pi,fo) ■= D 2 X1 + Dl 2 + (D X3 + p 2Xl - (3 lX2 )\ (3.3) 

in L 2 (R 3 + ). We define 

a{6) :=inf Spec P(sin6, cos 9). (3.4) 

We will need a number of properties of this spectral function. For proofs and more 
precise results see |HeMo3i ILuPa4j . 
Proposition 3.1. 

The function a(9) defined on [— 7r/2,7r/2] is continuous and even and satisfies 
• a is strictly increasing on [0, 7r/2]. 
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• o-(O) = 6 and <t(it/2) = 1. 

• We have the following expansion for small 9: 

o-(9) = Q + ai \e\+O(6 2 ), (3.5) 
where a% — v^o > (with So from fl3.2|) ). 

3.1.3. Montgomery's model. 

We finally mention briefly a model on L 2 (M. 2 ) which has been studied in |Mon[ 
IHeMoll IPaKwl IHe2| . This model has magnetic field vanishing on a line. Define 
Pm to be the self-adjoint realization of 

P M := Dl + (D y - x 2 ) 2 , 

on L 2 (R 2 ). Define the constant v$ as 

vq := inf SpecPjvi- (3-6) 

3.2. Localization estimates. 

We start by recalling the decay in the direction normal to the boundary. We will 
often use the notation 

t(x) :=dist(x,9fi). (3.7) 

Now, if <j) £ Co°(f2), i.e. has support away from the boundary, a simple standard 
integration by parts (already used in a very similar way in (|2.19| ) implies that 

Qb(4>) > B\\m (3.8) 

It is a consequence of this elementary inequality (and the fact that 0o < 1) that 
ground states are exponentially localized near the boundary. We give the result 
without proof since this has been proven by many authors in this and similar 
situations. 

Theorem 3.2. 

Let fl C K 3 be a bounded open set with smooth boundary. Then there exist constants 
C, a\ > 0, Bq such that 

e 2a lB ^t( x ) ^ b{x) \2 + B -i|(_iv + BF)ip B (x)\ 2 ) dx < CU B \\l (3-9) 

for all B > Bq, and all ground states ips of the operator Ti(B). 

We will mainly use this localization result in the following form. 
Corollary 3.3. 

Let Q, C M 3 be a bounded open set with smooth boundary. Then for all n eN there 
exists C n > such that 

J t{x) n \^ B (x)\ 2 dx < C n B- n ' 2 U B \\l 

for all B > and all ground states ips of the operator TL{B). 

We now define tubular neighborhoods of the boundary as follows. For e > 0, 
define 

B{dn, e) := {x € ft : t(x) < e}. (3.10) 
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For sufficiently small eo we have that for all x G B(dVL, 2eo) exists a unique point 
y = y(x) £ dfi such that t(x) = dist(x,y). We fix such an eo in the rest of the 
paper. 

Define, for y S dtt, the function 8(y) S [— ir/2,ir/2] by 

sin %):=-/?• iV(y). (3.11) 

We extend 9 to the tubular neighborhood B(dVt, 2cq) by 6{x) := 6(y(x)). 

In order to obtain localization estimates in the (boundary) variable normal to 
r, we use the following operator inequality (see |HeMo5[ Theorem 4.3]). 

Theorem 3.4. 

Let Q C M 3 be a bounded open set with smooth boundary. Let Bq be chosen such 
that B Q 3 ^ 8 = eo and define, for B > Bq, C > and x € Q, 

„ r , , \B-CBV\ t(x)>2B- 3 / s , ,„ . 

WB ' C(X) := (fl^x)) - ^)<2S-3/B. (3 ' 12) 

Then, there exists Co smc/i i/iai 

> W Bt c, (3.13) 
f in i/ie sense of quadratic forms ) for all B > Bq and all C > Cq . 

We use Theorem 13.41 to prove Agmon type estimates on the boundary. 
Theorem 3.5. 

Suppose that il C K 3 satisfies Assumptions [7751 |J.7| and \1.8l Define for x € dCl, 

d r (x) := distgn(a;,r), 

where distgn is the distance associated to the induced metric of d£l, and extend dr 
to a tubular neighborhood of the boundary by dr(x) := dr(y(x)), where y{x) is the 
unique boundary point closest to x. 
Then there exist constants C, 0,2 > 0,Bq, such that 

e ^B^ M ^ lMx)l 2 dx < Ce CBl "\\Ml (3-14) 

B(dn,e ) 

for all B > Bq and all ground states tpB ofH(B). 

It is useful to collect the following easy consequenc^l . 
Corollary 3.6. 

Suppose that ft C K 3 satisfies Assumptions [7751 U- "A and \1.8[ Then for all n 6 N 
there exists C n > such that 

d r (x) n \^ B (x)\ 2 dx < C n B- n ^U B \\l (3.15) 

B(an,e ) 

for all B > and all ground states t/jg ofH(B). 



° It is enough to consider separately the cases when dp < EB 4 and the case when dr > 
_ 1 

EB 4 with E large enough. 
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Proof of Theorem \3.5[ 

We may clearly assume that HV'slh = 1- Let xiiX2 £ C°°(R) such that Xi is 
decreasing, xi = on [2,+oc), xi = 1 on (— oo, 1], X2 is increasing, xi = 1 on 
[2,+oo), xi = on (-oo, 1]. 

By the standard localization formula we find, since TL{B)ipB = Xi(B)ipB, 

MB) X i{ L )X2{DB^d T )e^ Bl/2d r^ L ' 
eo 



= Qb \xi(-)X2(DB 1 ' 4 d T 
L e 



/ V(xi(-)x 2 (^ 1/4 dr) e a2Bl/2 ^ 2 ) 2 \M 2 dx. 
J e 



(3.16) 



Here D is a sufficiently small constant to be determined. 
We estimate using Theorem 13.41 

t 



Xl (l) X2 (DB^d r )e^ Bl/2 4 /2 ^ 



eo 



> 



/ W B {x) xi(-)X2pS 1/4 dr)e a2Bl/2d r / V s * dx. (3.17) 
J £Q 



Also, 



V(xi(-)X2(i? J B 1/4 dr) 



2 B 1/2 4 /2 N 



< 2 



Vxi(-) X^(^ 1/4 rfr)e 2a2Bl/2 ^ 2 
2|Vx 2 (^ 1/4 dr)| 2 X 1 (-)e 2a2Bl/2d r /2 



eo 



^a^d r x?(-)x^(^ 1/4 rfr)e 2 - Bl/2rf r /2 . 
Z e 



(3.18) 



Combining, (j3~TtJj) . (j3~T7) and l|3"35]) . we find 

/ (WB(aO - \i(B) - ^a^d r ( a; ))x 1 (-)x^^ 1/4 dr)e 2 - Bl/2d r /2 |^|2 da; 
J z e 



< C 



I B{dn,2e )\B(dn,e ) 
(7^ 1 / 2 e 4 v / 2 a 2- D ^ B * 



Z 2a * Bl/2d r 2 \4> B \ 2 dx 



{x£B{dn,2t ) ■ B 1 / i d r (x)<%} 



\ip B (x)\ 2 dx. (3.19) 



Since fl is bounded there exists Dx > such that dr(^) < Di for all X. Thus we 
can estimate, with a\ being the constant from Theorem [ 



B{dnat )\B(dil,co) 



< (7 e 2B 1/3 (a 2 0?/ 2 - ai£ o) ||^ B ||2 = 0( B -oo) s (3.20) 

3 

where the last estimate holds if ai is sufficiently small {a^D^ < aie$). 



ON THE GINZBURG-LANDAU CRITICAL FIELD IN THREE DIMENSIONS 



17 



Notice now that the assumption (|1.26j) implies that j3 ■ N vanishes exactly at 
order 1 on T. Therefore, using the boundedness of f2, there exists a constant C > 
such that 

C -1 d r (a?) < 0(x) < Cdr(x), (3.21) 
for all x G B(dft, 2eo). Therefore, we observe that 

• a is monotone, 

• ol\ > (from (|3~5")0 , 

• D can be chosen arbitrarily small, 

• Ai(-B) is bounded from above (see (|1.27p ). 

and we find that, if a 2 is sufficiently small, there exists Bq such that 

W B (x) - Ax(B) - ~<%Bd r (x) > B 3/4 , (3.22) 

for all x G B(dQ, 2e ) with d v {x) > D~ 1 B~ 1 / 4 , and B > B . 

Inserting (|3.20[) and (|3.22[) in (|3. 19|) . yields, for some constants C > and B , 

x \{L) x l {D B^U T )e 2a - Bl,2d r"\^ B ^dx < Ce CBi , VB > B . (3.23) 
£0 

Since e 2a * Bl/2d r 2 < e a 2 2^D-^ B i when B i/4 dr < ^ foUows from ([X231 . 

□ 

Consider now the set A^r C T where the function 70 is minimized, 

Mr := {x G T : 70 = 70}- (3.24) 

Theorem 3.7. 

Let (5 > and suppose that Q c R 3 satisfies Assumptions [7751 |j. 7| and[77JJ Then, 
for all N > 0, i/iere exists Cn > 0, swc/i i/iai if ipB is 1 normalized ground state 
eigenfunction ofTi(B), then 

f \iP B (x)\ 2 dx<C N B~ N , (3.25) 

/or aZZ B > 0. 

The proof of the theorem is based on a careful reading of what is proved, in 
[HeMo51 Section 150 in the second zone T which corresponds to the neighborhood 
of T defined by 

T 11 :={xen\ d r (x) + d(x, dft) < B~ s } 

whert0 S G]j|, §[. 

In this region it is shown that we have, for some rj > 0, the inequality 

B* I (70(2;) -%)\ijj B (x)\ 2 dx 
Jr 11 

< {{H{B)-\{B))^ B \^ B )+CB^U B \\ 2 . (3.26) 

Then we can proceed like in the previous proof. The only point here is that the 
situation is easier because we need only to have an estimate in the region of T 11 
where dist(x, .Mr) > <5 > (with 5 independent of B). 

4 with the correspondence h = B~ 1 , r(h) = h s . 
5 Note that | < jg. 
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It is clear that, outside this zone, we have already shown that ipB is exponentially 
small. 

Proposition 3.8. 

Let dr be the function defined in Theorem \3.5\ and let Tj be one of the curves making 
up T. Let sq € Tj and define, for e > 0, 

n(e, s ) = {x S ft : dist(x, T) < e and dist(x, s ) > e}. (3.27) 

Then, if e is sufficiently small, there exists a function <j) € C°°(f2) such that the 
magnetic potential A := F + V(j>, satisfies 

\A(x)\<c(t(x) + d r (x) 2 y 

for all x € f2(e, So). 

An easy localization argument shows that we can carry out the above gauge 
change simultaneously at each Tj. 

Corollary 3.9. 

Let (si, . . . , sjv) 6 Ti x ■ • • x IV and define, for e > 0, 

0(e, (si, . . . , sjv)) = {a; £ 51 : dist(a;, r) < e and mindist(x, Sj) > e}. (3.28) 

Then, if e is sufficiently small, there exists a function <f> € C°°(f2) such that A := 
F + V0, satisfies 

\A(x)\ <c{t{x)+dv{xf), 
for all x E fife, (si, . . . , Sjy)) . 

Proof o/ Prop. [M 

We use the adapted coordinates (r, s,i) near Ij defined in [HeMo51 Chapter 8]. We 
briefly recall their properties. Let Tj be parametrized by arc-length as 

ir -I 

^S, 1 3 s^TJs) Gdn. 

ZTT 

Given a point x G Q sufficiently close to Tj there exists a unique point y(x) £ 
dQ such that dist(x, dQ) — dist(a;, y(x)) and a unique point Tj{s{x)) € Tj such 
that dr{y{x)) — distgfi(y(x), Tj) = diston(y(x), Tj(s(x))). The coordinates (r, s,t) 
associated to the point x now satisfy 

|r| = dr(y[x)), s = s(x), t = dist(x, dfl). 

Let A\dr + A 2 ds + A^dt be the magnetic one-form uja = A • dx pulled-back 
(or pushed forward) to the new coordinates (r,s,t). Also write the corresponding 
magnetic two-form, duiA, as 

Bi 2 dr A ds + B^dr A dt + B 2 zds A dt. 

Clearly, 

B i:i = diAj - djli, (3.29) 
for i < j and where we identify (1, 2, 3) with (r, s, t) for the derivatives. 
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The magnetic field [3 corresponds to the magnetic two- form via the Hodge- map. 
In particular, since (3 is tangent to dQ. at L we get that 

5 12 (0,s,0) =0. (3.30) 
We now find a particular solution A to the equations (|3.29j) . We make the Ansatz 

Ai = - [ B 13 (r,s,T)dT, (3.31) 
Jo 

A 2 = - [ B 23 (r,s,T)dT + Mr,s), (3.32) 
Jo 

A 3 = 0. (3.33) 

Using the relation d{duj^_) = 0, we see that the above Ansatz gives a solution if ip 2 
is chosen as 



Mr,s)= / B 12 (p,s,0)dp. (3.34) 
Jo 

One verifies by inspection that with these choices 

\A\<C(r 2 +t). (3.35) 

By transporting this A back to the original coordinates we get the existence of an 
A with 

curlA=l, |A»| < C(t(x) +d r (x) 2 ). 

Since f2(e, sq) is simply connected (for sufficiently small e) A is gauge equivalent to 

F and the proposition is proved. □ 



3.3. The derivative. 

In this section we prove how one can derive the monotonicity result from the known 
asymptotics of the ground state energy and localization estimates for the ground 
state wave function itself. 

Based on these estimates the proof of Theorem 11.111 is very similar to the two- 
dimensional case. 
Proof of Theorem 

Let Ffi be the vector potential defined in (|1.15j) . In this proof we will abuse notation 
slightly and use the symbol 7~t(B) for the operator P%p n with Neumann boundary 
conditions. This corresponds to a gauge change (i.e. a unitary transformation) 
with respect to the previous notation. In this notation we have 

V{H{B)) = {ue W 2 > 2 {n) : N • Vu| ao = 0}, 

in particular the domain is independent of B. Applying analytic perturbation 
theory to H(B) we get the existence of ±(B). 

Let L = UjL^j be the decomposition of T in disjoint closed curves and let 
sj e Tj be a point with 'y(sj) > 70- Let f2(e, (s±, . . . , sjv)) be as defined in (|3 . 28[) 
with e so small that 

7(x) > 70 

for all x G 0(e, (si, . . . , sn)) ■ Let A be the vector potential defined in Corollarv l3.9l 
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Let Qb the quadratic form 

W 1,2 (Q) 3bh Qb{u) = [ \(-iVu+ BA)u\ 2 dx , 

and H(B) be the associated operator. Then Ti.(B) and Ti.(B) are unitarily equiva- 
lent: H(B) = e lB ' l 'H{B)e- lB < t ', for some </> independent of B. With ipf(-;/3) being 
a suitable choice of normalized ground state cigcnfunction, we get (by analytic 
perturbation theory applied to H(B) and the explicit relation between Tt(B) and 
H(B)), 

X' h+ (B) = (A>+( -,B), P B ^ti ■■:B)) + (p B ^t( ■ ;B) , Aipi( ■ ; B)) . (3.36) 
We now obtain for any b > 0, 

x[ A b) = Qb + ^U-,b))-q b (^(-,b)) _ b f |1|2 B)? dx 



n 



'^ -b \A\ 2 \^(^B)\ 2 dx. (3.37) 



Ai(S + 6)_-Ai(S) 
b 

We choose b := si -1 ^ with r/ < 77 (77 from ([L~2"7) ). Then, using (fL~2"7|) , ([337)1 

becomes 

V li+(S ) = 8 0+ ^-V»<l±^!^- CB -3 

-6 / |A| 2 (x; B)\ 2 dx , (3.38) 

for some constant C independent of B > Bo. If we can prove the existence of C 
such that 

B% [ |A| 2 \ifj+(x;B)\ 2 dx < C, (3.39) 

JQ. 

then we can take the limit B — > 00 in (|3.38j) and obtain 

liminf A' 1+ (B) > 6 . (3.40) 

B^oo ' 

Applying the same argument to the derivative from the left, _{B), we get 
(the inequality gets turned since b < 0) 

limsupAi _{B) < 6 . (3.41) 

Since, by perturbation theory, X[ + (B) < A^ _(£?) for all B, we get (jl.32p . 
Thus it remains to prove (|3.39|) . By Corollary 13. 91 we can estimate 

|A| 2 \rP+(x; B)\ 2 dx < C f (t 2 + r^+ix; B)\ 2 dx 

+ f , \i>t(x;B)\ 2 dx. (3.42) 

Jn\n(e,(s 1 ,...,s JV )J 



Combining Corollaries 13. 31 and 13. 61 and Theorem l3.7i we therefore find the existence 
of a constant C > such that : 

/ |A| 2 IV^x;^)! 2 ^ < CB- 1 , (3.43) 
Jn 

which is stronger than the estimate (|3.39[) needed. □ 
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Remark 3.10. 

Actually, since (I3.43|) is stronger than (|3.39|) . we do not need the full asymptotics 
(|1.27p but only the weaker result 

Ai(.B) = e Q B + o(B). 

Using this information only and choosing b = r]B instead (and taking the limit 
rj — ► in the end) we can still achieve the proof of Theorem \l.lli However, the 
second term in the asymptotics (|1.27[) is used in the proof of Theorem \3.7\ so there 
seems no gain in not using it at this point. 
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